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Abstract 



We apply a quasi-model-independent strategy ("Sleuth") to search for new high pr physics in w 
100 pb"'^ of pp collisions at -/s = 1.8 TeV collected by the D0 experiment during 1992-1996 at the 

Fcrmilab Tovatron. Over thirty-two e/iX, l^+jcts-likc, Z+jcts-likc, and {l/^)(l/'y)(l/'^)X exclusive 
final states are systematically analyzed for hints of physics beyond the standard model. Simultaneous 
sensitivity to a variety of models predicting new phenomena at the electroweak scale is demonstrated 
by testing the method on a particular signature in each set of final states. No evidence of new high pr 
pliA'sics is observed in the course of this search, and we find that 89% of an ensemble of hypothetical 
similar experimental runs would have produced a final state with a candidate signal more interesting 
than the most interesting observed in these data. 
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I. INTRODUCTION 

The standard model is an impressive theory, accurate- 
ly predicting, or at least accommodating, the results of 
nearly all particle physics experiments to date. It is gen- 
erally accepted, however, that there is good reason to 
believe that hints of new physics are likely to appear at 
or around the energy scale of 1 TeV. 

Electroweak symmetry is broken in the standard model 
when a scalar field (the Higgs field) acquires a vacuum 
expectation value. Since the quantum corrections to the 
renormalized mass squared of a scalar field grow as the 
square of the heaviest energy scale in the theory (naively 
the Planck scale, of order 10^^ GeV), and since the mass 
of the standard model Higgs boson is of the order of a 
few hundred GeV, a fine-tuning at the level of one part 
in 10^^ appears to be required to keep the Higgs mass at 
the electroweak scale. 

Two of the most popular solutions to this hierarchy 
problem are supersymmetry |]l| and strong dynamics [||. 
In their most general form these classes of models are 
capable of "predicting" any of many different signatures, 
depending upon the values that are chosen for the mod- 
el's parameters. Previous searches for these signals have 
fought to strike a balance between the simultaneous de- 
sires to assume as little as possible about the signal and 
yet achieve "optimal sensitivity" to more specific signals. 
These are necessarily contradictory objectives. 

Many new phenomena have been predicted in addition 
to those resulting from these proposed solutions to the hi- 
erarchy problem. Among them are leptoquarks, proposed 
in an attempt to explain the relationship between quarks 
and leptons in the standard model and appearing in many 
grand unified theories; composite quarks and leptons, in 
case the "fundamental" particles of the standard model 
turn out not to be fundamental at scales <, 10"^* meters; 
a fourth generation of quarks or leptons; excited quarks 
and leptons, in analogy to the excited states of hadrons 
observed at much lower energies; new heavy gauge bo- 
sons, arising from additional gauge symmetries in models 
extending the SU(3)cxSU(2)LxU(l)y of the standard 
model; and many others. Of course, Nature may have 
other ideas. The CDF and D0 collaborations have per- 
formed many searches on the data collected during Run 
I of the Fermilab Tevatron, but have we looked in all the 
right places? 

Figure ^ diagrams the final states that are populat- 
ed (i.e., that contain events) in the D0 Run I data. In 
this article we undertake a systematic and quasi-model- 
independent analysis of many of these exclusive final 
states, in the hope of finding some evidence for physics 
beyond the standard model. 
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FIG. 1. A diagram showing the final states populated in 
D0 data in Run I. Each row in a given column represents the 
final state defined by the objects in that row; to reduce clut- 
ter, jets are represented by an empty rectangle, rather than 
by a rectangle containing "j." Reading down the left col- 
umn are the final states e^lf^r, efifJrj, efifJr 2j, efifJr 3j, W, 
Wj, W2j, and so on. Rows with triangles (e.g., W and Wj) 
indicate final states analyzed previously by D0 in a man- 
ner similar to the strategy we use here, but without using 
Sleuth; rows with filled circles indicate final states analyzed 
with Sleuth. The remaining rows show populated final states 
not discussed in this article. 



tified with the D0 detector on an event-by-event basis, 
but taus are not; £ and the word "lepton" will therefore 
denote an electron (e) or a muon in this article. We 
use the composite symbol {£/j) to denote an electron, 
muon, or photon. X will denote zero or more objects, 
and (nj) will denote zero or more jets. Any inclusive fi- 
nal state [i.e., any state whose label includes the symbol 
X or (nj)] will refer to the physics objects actually re- 
constructed in the detector. Thus ee2j{nj) denotes the 
set of all events with two electrons and two or more jets. 
Any exclusive final state is defined according to the rules 
in Appendix |^. For example, since these rules include a 
prescription for identifying a Z boson from two charged 
leptons of the same flavor, we use ee 2j to denote the set 
of all events with two electrons and two jets having rriee 
substantially different from Mz, while events with two 
electrons and two jets having « Mz fall within the 
final state Z 2j. 

We begin in Sec. || by providing a brief review of the 
Sleuth search strategy and algorithm, and describing a 
slight change from the method advanced in Ref. [^. In 
Sec. Ill we discuss eight final states already analyzed by 



D0 in a manner similar to Sleuth, and motivate the fi- 
nal states to be considered in this article. In Sec. IV we 
describe the analysis of the W^+jets-like final states — 
events containing a single lepton, missing transverse en- 
ergy (^t), and two or more jets. In Sec. ^ we present 
the analysis of the Z-|-jets-like final states — events con- 
taining two leptons and two or more jets. In Sec. 



VI 



analyze the final states containing several objects, at least 
three of which are either a n el ectron, muon, or photon 
In Sec. VII we present the c ombined 
results of all of these final states. Section VIII contains 
our conclusions. 



II. Sleuth 



In Refs. we introduced a quasi-modcl-independent 
search strategy ("Sleuth"), designed to systematically 
search for new high px physics at any collider experi- 
ment sensitive to physics at the electroweak scale, and 
applied it to all events in the D0 data containing one or 
more electrons and one or more muons (efiX). Consid- 
ering again Fig. |l|, we see that the number of final states 
within efiX is a small fraction of the total number of fi- 
nal states populated by the D0 Run I data. If there is 
indeed a signal in the data, our chances of finding it grow 
proportionally to the number of final states considered. 

In this article we present a systematic analysis of 
thirty-two of these final states — those marked with a 
filled circle in Fig. ^ A large number of unpopulated fi- 
nal states with additional objects are analyzed implicitly; 
e.g., ee^lpT and efxflTj are among a host of unpopulated 
final states analyzed within the context of efiX . 

The notation we use to label final states may require 
explanation. Electrons and muons are confidently iden- 



In this section we provide for completeness a brief 
overview of the Sleuth algorithm, which is described in 
detail in Ref. and its application to the final states 
e^X. 

A. Search strategy 

We partition our data into exclusive final states, using 
standard identification criteria to identify electrons, mu- 
ons, photons, jets, missing transverse energy, and W and 
Z bosons. Although experimental realities will occasion- 
ally force slight modifications to these criteria, a set of 
standard definitions determined a priori is used wherever 
possible. 

The production and subsequent decay of massive, non- 
standard-model particles typically results in events con- 
taining objects with large transverse momentum (j>t)- 
For each exclusive final state we therefore consider the 
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small set of variables defined by Table |. In order 
to reduce backgrounds from QCD processes that pro- 
duce extra jets from gluon radiation, or two energetic 
jets through a t-channel exchange diagram, the notation 
^'py is shorthand for p!^ if the final state contains only 
one jet, 'Y^^=2Pt final state contains n > 2 jets, 

and X]r=3 Pt final state contains n jets and noth- 

ing else, with n > 3. Leptons and missing transverse 
energy that are reconstructed as decay products of W or 
Z bosons arc not considered separately in the left-hand 
column. Thus the variables corresponding to the final 
state Wjj, for example, are and X)'pt' Pt ^^^^ 
are not used, even though the events necessarily contain 
a lepton and missing transverse energy, since the lep- 
ton and missing transverse energy have been combined 
into the W boson. Since D0's muon momentum reso- 
lution in Run I was modest, we define X^-Pt ~ TIPt 
for events with one or more electrons and one or more 
muons, and we determine the missing transverse energy 
from the transverse energy summed in the calorimeter, 
which includes the pT of electrons, but only a negligible 
fraction of the pt of muons. When there are exactly two 
objects in an event (e.g., one Z boson and one jet), their 
Pt values are expected to be nearly equal, and we there- 
fore use the average pr of the two objects. When there 
is only one object in an event (e.g., a single W boson), 
we use no variables, and simply perform a counting ex- 
periment. We expect evidence for new physics to appear 
in the high tails of these distributions. 



If the final state includes 



then consider the variable 



one or more charged leptons 
one or more electroweak bosons 
one or more jets 



TABLE I. A quasi-model-independently motivated list of 
interesting variables for any final state. The set of variables 
to consider for any particular final state is the union of the 
variables in the second column for each row that pertains to 
that final state. 



B. Algorithm 

Although the details of the algorithm are complicated, 
the concept is straightforward. What is needed is a data 
sample, a set of events modeling each background process 
i, and the number of background events bi ± Sbi from 
each background process expected in the data sample. 
From these we determine the region of greatest excess and 
quantify the degree to which that excess is interesting. 

The algorithm, applied to each individual final state, 
consists of seven steps. 

1. We begin by constructing a mapping from the d- 
dimcnsional variable space defined by Table Q into 



the d-dimensional unit box (i.e., [0, 1]'') that flat- 
tens the background distribution, and we use this 
to map the data into the unit box. This change 
of variable space greatly simplifies the subsequent 
analysis. 

2. Central to this algorithm is the notion of a "region" 
about a set of 1 < < A^data data points, defined 
as the volume within the unit box closer to one of 
the data points in the set than to any of the other 
data points in the sample. The arrangement of data 
points themselves thus determines the regions. A 
region containing A^ data points is called an A^- 
region. 

3. Each region R contains an expected number of 
background events bn, equal to the volume of the 
region x the total number of background events 
expected, and an associated systematic error Sb^, 
which varies within the unit box according to the 
systematic errors assigned to each contribution to 
the background estimate. We can therefore com- 
pute the probability p^ that the background in 
the region fluctuates up to or beyond the observed 
number of events. This probability is our first mea- 
sure of the degree of interest of a particular region. 

4. The rigorous definition of regions reduces the num- 
ber of candidate regions from infinity to « 2^''''''' . 
Imposing explicit criteria on the regions that the al- 
gorithm is allowed to consider further reduces the 
number of candidate regions. (See Sec. IIP .) Our 
assumption that new physics is most likely to ap- 
pear at high Pt translates to a preference for re- 
gions in a particular corner of the unit box; criteria 
are thus constructed to define "reasonable" discov- 
ery regions. The number of remaining candidate 
regions is still sufficiently large that an exhaustive 
search is impractical, and a heuristic is employed 
to search for regions of excess. In the course of this 
search the A^-region TZ n for which p^ is minimum 
is determined for each A^, and pn = min^j (p^) is 
noted. 

5. In any reasonably-sized data set, there will always 
be regions in which the probability for bn to fluctu- 
ate up to or above the observed number of events 
is small. The relevant issue is how often this will 
happen in an ensemble of hypothetical similar ex- 
periments (hse's). This question can be answered 
by performing these hse's; i.e., generating random 
events drawn from the background distribution, 
and computing pN by following steps (l)-(4). Gen- 
erating many such hse's, we can determine the frac- 
tion P/v of hse's in which the pn found for the hse 
is smaller than the pN observed in the data. 

6. We define P and A^min by P = P/v„i„ = minjv (Pn), 
and identify TZ = ^Ar„i„ as the most interesting 
region in this final state. 
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Data set 


Fakes 


Z ^ TT 


7* TT 


WW 


tt 


Total 


Data 




18.4±1.4 


25.6±6.5 


0.5±0.2 


3.9±1.0 


0.011±0.003 


48.5±7.6 


39 




8.7±1.0 


3.0 ±0.8 


0.1±0.03 


1.1±0.3 


0.4±0.1 


13.2±1.5 


13 




2.7±0.6 


0.5±0.2 


0.012±0.006 


0.18±0.05 


1.8±0.5 


5.2±0.8 


5 




0.4±0.2 


0.07±0.05 


0.005±0.004 


0.032±0.009 


0.7±0.2 


1.3±0.3 


1 




30.2±1.8 


29.2±4.5 


0.7±0.1 


5.2±0.8 


3.1±0.5 


68.3±5.7 


58 



TABLE II. The numbers of expected background events for the populated final states within e/iX. The uncertainties in e\iX 
are smaller than in the sum of the individual background contributions obtained from Monte Carlo because of an uncertainty 
in the numbers of extra jets arising from initial and final state radiation in the exclusive channels. 



7. We use a second ensemble of hse's to determine the 
fraction V of hse's in which P found in the hse is 
smaller than P observed in the data. The most im- 
portant output of the algorithm is this single num- 
ber V , which may loosely be said to be the "fraction 
of hypothetical similar experiments in which you 
would see something as interesting as what you ac- 
tually saw in the data." V takes on values between 
zero and one, with values close to zero indicating 
a possible hint of new physics. In computing V 
we have rigorously taken into account the many re- 
gions that have been considered within this final 
state. 

The smallest V found in the many different final states 
considered (T^min) determines "P, the "fraction of hy- 
pothetical similar experimental runs (hser's) that would 
have produced an excess as interesting as actually ob- 
served in the data," where an hser consists of one hse for 
each final state considered. V is calculated by simulating 
an ensemble of hypothetical similar experimental runs, 
and noting the fraction of these hser's in which the small- 
est V found is smaller than Vmin- The correspondence 
between V and Vmin is determined to zeroth order by the 
number of final states considered in which the expected 
number of background events is ^ 1, with "smaller" final 
states contributing first order corrections. V also takes 
on values between zero and one, and the potential pres- 
ence of new high px physics would be indicated by finding 
V to be small. The difference between V and V is that in 
computing V we account for the many final states that 
have been considered. V can be translated into units of 
standard deviations [Via]) by solving the unit conversion 
equation 



pendix B 1 contains examples of the types of new physics 



V 



1 



2tt Jp 



(1) 



for Vta]- A similar equation relates V and Vta]- 



C. efj,X 

In Ref. 1^ we applied Sleuth to the e^X final states, 
using a data set corresponding to 108±6 pb"^ of inte- 
grated luminosity. We summarize those results here. Ap- 



that might be expected to appear in these final states. 

Events containing one or more isolated electrons and 
one or more isolated muons, each with px > 15 GeV, 
are selected. Global cleanup cuts are applied to remove 
events in which there was activity in the Main Ring, the 
accelerator that feeds the Tevatron, reducing the total 
number of events by 30%. The dominant standard model 
and instrumental backgrounds to this data set are: 

• top quark pair production with t — > Wb, and with 
both W bosons decaying leptonically, one to ei^ (or 
to TV — > evw) and one to fj,v (or to tv — > ^vvv); 

• W boson pair production with both W bosons de- 
caying leptonically, one to ev (or to tv evuv) 
and one to ^.v (or to tv ^ivuv)] 



zh* 



e[ivvvv\ 



and 



• instrumental ( "fakes" ) : W production with the W 
boson decaying to [iv and a radiated jet or photon 
being mistaken for an electron, or bb/cc production 
with one heavy quark producing an isolated muon 
and the other being mistaken for an electron ^ . 

The numbers of events expected for the various samples 
and data sets in the populated final states within e^X 
are given m Table ||. 

Among the systematic errors in these and other final 
states is an uncertainty in the modeling of additional ra- 
diated jets. Our consideration of exclusive final states 
makes this error more important than if inclusive final 
states were considered. An uncertainty of « 20% in the 
number of expected events, obtained by comparing the 
jets radiated by various Monte Carlo programs, is added 
in quadrature to systematic errors from other sources to 
obtain the total systematic error quoted in Table |^ and 
elsewhere. Because final states are analyzed indepen- 
dently, and because the definition of V depends only on 
the smallest V found, we can, to first order, ignore the 
correlations of uncertainties among different final states. 

We demonstrated Sleuth's sensitivity to new physics 
by showing that the method is able to find indications 
of the existence of WW and ti production in these final 
states when the backgrounds are taken to include only 
Z/7* TT and fakes. Figure || shows our sensitivity to 
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tt in an ensemble of mock data samples when the back- 
grounds include WW in addition to Z/^* — > tt and 
fakes. All samples with V[cr] > 2.0 appear in the right- 
most bin. We see that Sleuth, with no knowledge of the 
top quark's existence or characteristics, finds Via] > 2.0 
in over 25% of the mock samples. (For mock samples 
containing only Z/j* tt, fakes, and WW, the dis- 
tribution is roughly Gaussian and centered at zero with 
unit width.) After performing these sensitivity checks, we 
added all known standard model processes to the back- 
ground estimate and searched for evidence of new high 
Pt physics. The result of this analysis is summarized in 
Table III , No evidence of new physics is observed. 



c 

cr 
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Backgrounds: WW, Fakes, Z/y'—>TT 
Mock Samples; tt, WW, Fakes, Z/y'— >tt 

■ AntiCornerSphere+ Isolation 

Hyperplanes+Connectivity+ReasonableSize 

Combination of e/i^r. s/j.^^], e^i^, 2J and e/i^r 3j 
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-1 
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FIG. 2. Distribution of Vi^r] in an ensemble of mock exper- 
imental runs on the four exclusive final states e/i-^r, e^lprj, 
efi^T 2j, and e/i-^r 3j. The background includes Z/y* — > tt, 
fakes, and WW. The mock samples making up the distribu- 
tions contain tt in addition to Z/"/* — > tt, fakes, and WW. 



D. Region criteria 

Use of Sleuth requires the specification of criteria that 
define the regions that Sleuth is allowed to consider. In 
the analysis of e^X we imposed two criteria: AntiCorner- 
Sphere (c^), which restricts the allowed region to be de- 
fined by those data points greater than a distance r from 
the origin of the unit box, where r is allowed to vary; 
and Isolation (c^), which requires that there exist no da- 
ta points outside the region that are closer than ^ to any 
data point inside the region, where ^ = ^/{^N^itg) is a 
characteristic distance between the A'data data points in 
the d-dimensional unit box. 

For the analysis described in this article we use Hyper- 
planes (c^), a criterion defined but not used in Ref. 
Hyperplanes is less restrictive than AntiCornerSphcre, 
in the sense that any region satisfying AntiCornerSphcre 
will also satisfy Hyperplanes. Hyperplanes has the ad- 
vantage of allowing regions that lie in the high tails of 
only a subset of the variables considered. A region R 
in a d-dimensional unit box is said to satisfy Hyper- 
planes if, for each data point p inside R, one can draw a 
{d — l)-dimensional hyperplane through p such that all 
data points on the side of the hyperplane containing the 
point 1 (the "upper right-hand corner of the unit box" ) 
are inside R. An example of a region satisfying Hyper- 
planes is shown in Fig. 

We continue this boolean criterion to the unit inter- 
val [0, 1] in order to ensure the continuity of the final 
result under small changes in the background estimate. 
For each data point i inside the candidate region R and 
each hyperplane hi through i, we define djhi to be the 
distance between a data point j lying outside R and the 
hyperplane hi. This quantity is taken to be positive if j 
and the point 1 are on the same side of hi, and negative 
otherwise. Letting 



9{x) 



,x <0 

X ,Q<x <l 

1 , 1 < a; 



(2) 



we define 



Data sot 


V 




0.14 




0.45 




0.31 


ep,lpT 3j 


0.71 



TABLE III. Summary of results on all final states within 
ep,X when all standard model backgrounds, including tt, are 
included. We note that all final states within e^X have been 
analyzed, including (for example) ee^Ipr and enfrj. All 
final states within efiX but not listed here are unpopulated, 
and have P = 1.00. 



min max djh- /£,). 

hi j^R 



(3) 



Loosely speaking, the introduction of c|| corresponds to 
widening the lines drawn in Fig. ^ into bands of width ^, 
choosing c|[ = 1 if all data points "up and to the right" of 
these bands are inside R, finding = if there is a point 
"up and to the right" that is not inside R, and choosing 
c|j between and I if there are one or more points not 
inside R lying on the bands. Note that reduces to the 
boolean operator of the preceding paragraph in the limit 
^ 0, corresponding to the squeezing of the bands back 
into the lines in Fig. ^. 



We also impose the criterion Connectivity 



to en- 



sure connected regions, and the criterion ReasonableSize 
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Symbol 
— — 



Name 



A region satisfies this criterion if 




FIG. 3. An example of a region satisfying Hyperplanes. 
The boundary of the figure is the unit box; open squares repre- 
sent data points outside the region R\ filled squares represent 
data points inside the region R. The three dashed lines indi- 
cate hyperplanes hi (which are lines in this two-dimensional 
case) that can be drawn through the points at (x, y)i — 
(0.34,0.96), (0.74,0.95), and (0.935,0.515) with the proper- 
ty that all of the data points "up and to the right" of hi are 
inside R. 



AntiCornerSphere One can draw a sphere centered 
on the origin of the unit box con- 
taining all data events outside the 
region and no data events inside 
the region. 

(/ Isolation There exist no data points outside 

the region that are closer than ^ to 
any data point inside the region. 
Hyperplanes For each data point p inside R, 
one can draw a {d— l)-dimensional 
hyperplane through p such that 
all data points on the side of the 
hyperplane containing the point 1 
are inside R. 

Connectivity Given any two points a and b with- 
in the region, there exists a list of 
points pi = a,p2, • . • ,p„-i,p„ = 6 
such that all the Pi are in the re- 
gion and Pi+i is a neighbor of pi. 
ReasonableSize The region contains fewer than 50 
data points. 

TABLE IV. Summary of the region criteria imposed in our 
previous analysis of efiX (above middle line) and those im- 
posed in the analyses described in this article (below middle 
line). 1^ = l/(4A'^^^t''^) is a characteristic distance between the 
A'^data data points in the d-dimensional unit box. 



highest energy scales currently achievable, allows for the 
possibility that something there may yet remain undis- 
covered. 



(c^) to limit the size of the regions we consider to that 
expected for a typical signal and to reduce the compu- 
tational cost of finding the most interesting region. A 
region R is said to satisfy Connectivity if, given any two 
points a and b within R, there exists a list of points 
Pi = a,p2, . ■ . ,Pn-i,Pn = b such that all the pi are in R, 
and the 1-region about Pi+i shares a border with the 
1-region about pi. A region is said to satisfy Reason- 
ableSize if it contains fewer than 50 data points. These 
criteria are summarized in Table IV . 

In Ref. [H we demonstrated Sleuth's ability to find in- 
dications of ti in the e^X final states using the criteria 



Figure H shows that the combination c 



(sol- 



id) performs similarly to those criteria (dashed) in this 
test. 



III. CHARTED AND UNCHARTED TERRITORY 

The D0 experiment |^ began collecting data at ^/s ~ 
1.8 TeV in 1992, and completed its first series of runs in 
1996. These data have been carefully scrutinized by the 
D0 Collaboration. Nonetheless, the incredible richness 
of these data, which probe fundamental physics at the 



A. Final states already considered by D0 

Some portions of these data have been more compre- 
hensively scrutinized than others. In particular, there are 
eight final states — those marked with triangles in Fig. |l| 

— that D0 has already analyzed in a manner similar to 
the Sleuth prescription. 

In final states that contain only a single object (such 
as a W or Z boson), there are no non-trivial momen- 
tum variables to consider, and the Sleuth search strategy 
reduces in this case to a counting experiment. In final 
states containing exactly two objects (such as ee, Zj, 
or Wj), the single momentum variable available to us is 
the average (scalar) transverse momentum of the two ob- 
jects, assuming that both are sufficiently central. D0 has 
analyzed eight final states in these limiting cases. These 
analyses do not precisely follow the Sleuth prescription 

— they were performed before Sleuth was created — so 
V is not calculated for these final states. Nonetheless, 
they are sufficiently close to our prescription (and there- 
fore sufficiently quasi-model-independent) that we briefiy 
review them here, both for completeness and in order to 
motivate the final states that we treat in Sees. ^-[v|. 
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Examples of the types of new physics that could be ex- 
pected to ap pear in a few of these final states are provided 
in Appendix 



B2 



a. 2j. D0 has performed an analysis of the dijet 
mass spectrum and angular distribution |8| in a search 
for quark compositeness. We note that the dijet mass 
and the polar angle of the jet axis (in the center-of- 
mass frame of the system) together completely charac- 
terize these events, and that two central jets with large 
invariant mass also have large average px ■ No compelling 
evidence of an excess at large jet transverse momentum 
is seen in either case. 

b. W. The Sleuth-defined W final state contains all 
events with cither: one muon and no second charged lep- 
ton; or one electron, significant missing transverse energy, 
and transverse mass 30 < < 110 GeV. The Sleuth 
prescription reduces to a cross section measurement in 
this case. D0 has measured the inclusive W boson cross 
section , and finds it to be in good agreement with the 
standard model prediction. 

c. ef!T- Events that contain one electron, no second 
charged lepton, substantial fx, and have transverse mass 

> 110 GeV belong to the efr final state. This final 
state contains two objects (the electron and the missing 
transverse energy), so we consider the average object pr, 
which is approximately equal in this case to /2. D0 
has performed a search for right-handed W bosons and 
heavy W bosons in 79 pb~^ of data jl^], looking for an 
excess in the tail of the transverse mass distribution. No 
such excess is observed. 

d. Wj. In the two-object final state Wj, the average 
transverse momentum of the two objects is essentially 
p^ , the transverse momentum of the W boson. D0 has 
measured the W boson pT distribution |ll], and finds 
good agreement with the standard model. 

e. W^. Similarly, the transverse momentum distri- 
bution of the photon in WjX events has been analyzed 
by D0 in a measurement of the WW^ gauge boson cou- 
pling parameters |jl^. No excess at large p^ is observed. 
(The Sleuth prescription for defining final states is less 
well satisfied in D0's corresponding measurement of p^ 
in Z-fX events 

f. Z. As in the case of the W final state, our pre- 
scription reduces to a counting experiment in the Z final 
state. D0 has published a measurement of the inclu- 
sive Z boson cross section and finds it to be in good 
agreement with the standard model prediction. 

g. ee. Events containing two electrons and nothing 
else fall into the final state ee if the invariant mass rriee is 
outside the Z boson mass window of (82, 100) GeV. The 
single variable we consider in this two-object final state 
is the average scalar transverse momentum of the two 
electrons, which is simply related to the invariant mass 
iriee for sufficiently central electrons. D0 has analyzed 
the high mass Drell-Yan cross section in a search for indi- 
cations of quark-lepton compositeness with the full data 
set p^ , and has analyzed the ee invariant mass distri- 
bution in the context of a search for additional neutral 



gauge bosons in a subset of those data |T^. No dis- 
crepancy between the data and expected background is 
observed. 

h. Zj. In the two-object final state Zj, the average 
transverse momentum of the two objects is essentially the 
transverse momentum of the Z boson. D0's published 
measurement of the Z boson pT distribution p^ ] is in 
good agreement with the standard model prediction. 



B. Final states considered in this article 

The decision as to which of the remaining final states 
should be subjected to a Sleuth analysis was made on 
the basis of our ability to estimate the standard mod- 
el and instrumental backgrounds in each final state, and 
the extent to which a systematic analysis for new phys- 
ics is lacking in each final state. The final states we 
chose to analyze arranged themselves into four "classes" : 
efiX, VT-f jets-like final states, Z-t-jets-like final states, 
and {£ / / ■^){£ / j)X . The first of these clas ses h as been 
analyzed in Ref. |3| and summarized in Sec. II C 



A sys- 
asses of 



tematic Sleuth analysis of the remaining three c 
final states is the subject of the next three sections 



IV. JETS-LIKE FINAL STATES 

In this section we analyze the W^-|-jets-like final states 
— events containing a single lepton, missing transverse 
energy, and two or more jets. In Sec. [VA we describe 



the epT 2j{nj) and fi^ x '2j{n j) data sets and background 
estimates, and in Sec. IV B we present the results. Af- 



ter this, we feign ignorance of the heaviest quark in the 
standard model and check the sensitivity of our method 
to top quark pair production in Sec. IVC. A few of the 



many signals that mig ht a ppear in these final states are 
described in Appendix B3. 



A. Data sets and background estimates 

1. e$T2j{nj) 

The elflx '2.i{nj) data set [|l^ comprises 115±6 pb"^ of 
collider data, collected with triggers that require the pres- 
ence of an electromagnetic object, with or without jets 
and missing transverse energy. Offline event selection re- 
quires: one electron with transverse energy > 20 GeV 
and pseudorapidity |?7dct|< 1-1 or 1.5 <|T?dct|< 2.5 [ p^ ; 
IpT > 30 GeV; and two or more jets with pip > 20 GeV 
and |?7dct|< 2.5. Effects of jet energy mismeasurement 
are reduced by requiring the vector to be separated 
from the jets by > 0.25 radians if fr < 120 GeV. 
To reduce background from a class of events in which a 
fake electron's energy is overestimated, leading to spu- 
rious f^T, we reject events with pjf < 40 GeV. Events 
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Final State VK+jets QCD fakes ti Total Data 

e^T 2j 6.7 ± 1.4 3.3 ±0.9 1.7 ±0.6 11.6 ± 1.7 7~ 

e^rSj 1.0 ±0.4 0.48 ± 0.22 1.0 ± 0.4 2.5 ± 0.6 5 

e^T4j 0.15 ±0.11 0.38 ± 0.19 0.26 ±0.09 0.80 ± 0.24 2 

W{-^ e^r) 2j 334 ± 51 12.0 ± 2.6 4.0 ± 1.4 350 ± 51 387 

W{^e^T)Sj 57 ±9 3.4 ±0.9 6.0 ±2.1 66 ± 9 56 

W{^e^T)4,j 5.9 ±1.3 1.1 ± 0.4 3.9 ± 1.4 10.9 ±1.9 11 

W{~*e:^T)5j 0.8 ±0.3 0.19 ± 0.12 0.73 ±0.26 1.8 ±0.4 1 

W{-^e^T)Gj 0.12 ±0.06 0.030 ±0.015 0.10 ± 0.04 0.25 ±0.07 1 



TABLE V. Expected backgrounds to the ef^r 2j(nj) final states. The final states labeled "W(^ e^r)" have < 110 GeV; 
the final states labeled "e_^T" have > 110 GeV. We have extrapolated our background estimates to final states with five 
or more jets. Berends scaling and the data in this table suggest that a factor of ~ 7 in cross section is the price to be paid for 
an additional radiated jet with transverse energy above 20 GeV. 



Final State 


W+jets 


Z±jets 


WW 


tt 


Total 


Data 




48± 15 


1.6 ±0.4 


0.5 ±0.3 


0.42 ±0.14 


50± 15 


54 




10 ±3 


0.27 ±0.08 


0.41 ±0.26 


0.58 ±0.20 


11 ±3 


11 




2.8 ± 1.3 


0.022 ±0.011 




0.61 ±0.21 


3.5 ±1.3 


4 



TABLE VL Expected backgrounds for the P^t) 2j(nj) final states. 



containing isolated muons appear in a sample analyzed 
previously with this method (e/iX), and are not consid- 
ered here. 

The dominant standard model and instrumental back- 
grounds to the efr'^iinj) final states are from: 

• W -\- jets production, with W ev\ 

• multijet production, with mismeasured If^T and one 
jet faking an electron; and 

• tt pair production, with t — > Wb and with at least 
one W boson decaying to an electron or to a tau 
that in turn decays to an electron. 

The W^+jets background is simulated using VEC- 
BOS with HERWIG used for fragmenting the par- 
tons. The background from multijet events containing 
a jet that is misidentified as an electron, and with 
arising from the mismeasurement of jet energies, is mod- 
eled using multijet data. The probability for a jet to 
be misidentified as an electron is estimated to be 
(3.50 ± 0.35) X 10^''. The background from ti decays in- 
to an electron plus two or more jets is simulated using 
HERWIG with a top quark mass of 170 GeV. All Monte 
Carlo event samples are processed through the D0 de- 
tector simulation based on the GEANT [|2^ package. 

We estimate the number of events in the W^±jets-like 
final states to be 18 ±6 using the measured ti production 
cross section of 5.5 ± 1.8 pb j23|. The multijet back- 
ground is estimated to be 21 ± 7 events, using a sample 
of events with three or more jets with fx > 30 GeV. 
This is done by multiplying the fake probability by the 
number of ways the events satisfy the selection criteria 
with one of the jets passing the electron pT and require- 
ments. After the estimated numbers of ti and multijet 
background events are subtracted, the number of events 



with transverse mass of the electron and neutrino {m'^) 
below 110 GeV is used to obtain an absolute normaliza- 
tion for the VF±jets background. 

Following the Sleuth prescription, we combine the elec- 
tron and missing transverse energy into a W boson if 
30 < mf^ < 110 GeV, and reject events with < 
30 GeV. The expected numbers of background events for 
the exclusive final states within this efT'^jinj) sample 
are provided in Table 0. 

2. M^T 2j{nj) 

The ^^T^j^nj) data set [|4j corresponds to 94 ± 
5 pb""'^ of integrated luminosity. The initial sample is 
composed of events passing any of several muon ± jets 
triggers requiring a muon with pt^ > 5 GeV within 
|?7dot|< 1-7 and one or more jets with > 8 GeV and 
|77det|< 2.5. Using standard jet and muon identification 
criteria, we define a final sample containing one muon 
withpT > 25 GeV and |77dct|< 0.95, two or more jets with 
pip > 15 GeV and |77det|< 2.0 and with the most energet- 
ic jet within |77dot|< 1-5, and missing transverse energy 
-^T > 30 GeV. Because an energetic muon's momentum 
is not well measured in the detector, we are unable to 
separate 'W-like" events from "non-VF-like" events us- 
ing the transverse mass, as we have done above in the 
electron channel. The muon and missing transverse en- 
ergy are therefore always combined into a W boson. 

The dominant standard model and instrumental back- 
grounds to these final states are from: 

9 W + jets production with W 

• Z + jets production with Z —>■ fi^, where one of 
the muons is not detected; 
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• WW pair production with one W boson decaying 
to a muon or to a tau that in turn decays to a muon; 
and 

• tt pair production with t Wb and with at least 
one W boson decaying to a muon or to a tau that 
in turn decays to a muon. 

Samples oi W + jets and Z + jets events are gener- 
ated using VECBOS, employing herwig for parton frag- 
mentation. Background due to WW pair production is 
simulated with pythia Background from tt pair 

production is simulated using herwig with a top quark 
mass of 170 GeV. All Monte Carlo samples are again pro- 
cessed through a detector simulation program based on 
the GEANT package. 

The expected backgrounds for the exclusive final states 
within ^^T2j(nj) are listed in Table |v]. These W{-^ 
li^T)'^j{nj) final states are combined with th e W{—^ 
elpT) '^jinj) final states described in Sec. IV A 1| to form 

> 



Final State 



Total 



Data 



the W2j{nj) final states treated in Sec. IV A3 



consistency in this combination, we also require pj^ 
40 GeV for the W{-^ ^v) 2j{nj) final states. 



3. W2j{nj) 

Combining the results in Tables 0and0 gives the ex- 
pected b ackg rounds for the W2j{nj) final states shown 
in Table VII. We note the good agreement in all final 



states between the total number of background events 
expected and the number of data events observed. This 
of course is due in part to the method of normalizing 
the iy-|-jets background. The agreement in the final 
states containing additional jets is also quite good. A 
more detailed comparison between data and background 
in the more heavily populated final states {W2j, W3j, 
and W 4:j) is provided in Appendix 

Monte Carlo programs suitable for estimating back- 
grounds to final states with many additional jets are not 
readily available. It has been observed that the rate of 
a process may be related to the rate of the process with 
an additional radiated jet by a multiplicative factor of 
1/4-1/7, depending upon the pt and rj thresholds used 
to define a jet — this phenomenological law is known 
as Berends scaling |l9j. We estimate that this factor is 
« 1/5 for jets with |7;dct|< 2.5 andpr > 15 GeV, and that 
it is « 1/7 for jets with |77dct|< 2.5 and pr > 20 GeV. This 
will be used to estimate particular background contribu- 
tions to final states in which the expected background is 
<1. 



B. Results 

The results of applying Sleuth to the e^r 2j{n j) and 
W 2j(nj) data sets are summar ized in Table VIII and in 
Figs. ^ and |[ Recall from Sec. II B that the positions of 



W2j 
W3j 
WAj 
W5j 
W6j 



400 ± 53 
77± 10 
14.3 ± 2.3 
1.8 ±0.4 
0.25 ±0.07 



441 
67 
15 
1 
1 



TABLE VII. Expected backgrounds to the W 2j{nj) final 
states. 



Data set 



V 



e$T 2j 
e$T 3j 
elpT 4j 
W2j 
W3j 
WAj 
W5j 
W6j 



0.76 
0.17 
0.13 
0.29 
0.23 
0.53 
0.81 
0.22 



TABLE VIII. 
W2jinj). 



Summary of results on e^T 2j{nj) and 



the data points within the unit box are determined by 
the background distribution, which defines the transfor- 
mation from the original variable space, in addition to 
the location of the data points in that original space. We 
observe quite good agreement with the standard model 
in the VF±jets-like final states. 



C. Sensitivity checlt: tt ^+jets 

In this section we check Sleuth's sensitivity to in the 
final states W3j, W Aj, W5j, and W6j. After briefly 
putting this signal into context, we test Sleuth's ability 
to find tt in the data, and then in an ensemble of mock 
experiments. 

In 1997 D0 published a measurement of the top quark 
production cross section |Q based on events in the dilep- 
ton, ^±jets, ^±jets(//u), and "e^" channels, where "/a^" 
indicates that one or more of the jets contains a muon, 
and hence is likely to be the product of a 6 quark. 19 
events with no 6-quark tag are observed in ^±jets (nine 
events in the electron channel, and ten events in the muon 
channel) with an expected background of 8.7 ± 1.7. An 
additional eleven events are observed with a 6-quark tag 
(five events in the electron channel, and six events in the 
muon channel) with an expected background of 2.5 ± 0.5 
events. Three or more jets with px > 15 GeV are re- 
quired in both cases. The number of events observed in 
all four channels is 39 with an expected background of 
13 ± 2.2 events. The probability for 13 ± 2.2 to fluctuate 
up to or above 39 is 6 x 10"'', or 4.8 standard devia- 
tions. In the i+jets channel alone, the probability that 
8.7 ±1.7 fluctuates to the 19 events observed is 0.005, 
corresponding to a "significance" of 2.6a. 
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CL 




FIG. 4. The positions of the transformed data points in 
the final states elpr'^j, e^rSj, and e^T^j- The data points 
inside the region chosen by Sleuth are shown as filled cir- 
cles; those outside the region are shown as open circles. For 
these final states the variables py, ^t, and X/'Pt '"^^^ 
sidered, and the unit box is in this case a unit cube. The 
two-dimensional views shown here are the projections of that 
cube onto three orthogonal faces. 



Figures ^(a) and ^(c) show where tt Monte Carlo 
events fall in the unit box in the final states W 3j and 
W Aj. The distribution of these events is quite diffuse in 
the case of W3j, since tt is similar to the background in 
the variables pj^ and X^'pt ^^^^ channel. In the W Aj 
final state tt tends to populate regions of large pip, but 
the signal is nearly indistinguishable from background in 
the variable p}p . A check of Sleuth's ability to find tt in 
the W3j{nj) final states tests how well Sleuth performs 
when the signal shows up in a subset of the variables we 
choose to consider. 

Figures ^(b) and ^(d) show D0 data in the final states 
W 3j and WAj, when tt is not included in the back- 
ground estimate. Notice that the region chosen by Sleuth 
in the VF3j final state in Fig. ^(b) is very similar to 
the region populated by tt in Fig. ^(a). In the W 4j 
final state (d), the region chosen by Sleuth is nearly 
the entire unit box. Comparison with Fig. ^ shows 
how the absence of tt in the background estimate in 
this figure affects the transformation from the original 
variable space into the unit box. Applying Sleuth to 
these data while continuing to feign ignorance of ti, we 
find Vw3j = 0.12, Vw4j = 0.18, Vw5j = 0.37, and 
'PwGj — 0.09. Upon combining these results, we find 
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FIG. 5. The positions of the transformed data points in the 
final states W2j, W3j, W Aj, and Wbj. The data points 
inside the region chosen by Sleuth are shown as filled circles; 
those outside the region are shown as open circles. The single 
event in the W 5j final state is in the lower right-hand corner 
of the unit square, having X) ~ 300 GeV. 



■Pmin = min{Vw3j,'Pwij,'Pw5j,'Pw6j) = 0.09 (l.Sfj). 

Figure ^ shows a histogram of Vmin for a sample of 
mock experimental runs in which the backgrounds in- 
clude T4^-|-jets and QCD events, and the mock samples 
include tt in addition to the expected background. The 
number of background and tt events in the mock samples 
are allowed to vary according to statistical and systemat- 
ic errors. Note that since four final states are considered, 
the distribution of Vmin for an ensemble of experiments 
including background only has a median of « la. We 
see that Sleuth is able to find indications of the presence 
of ti in these final states, returning Vmin [a] > 3 in 30% 
of an ensemble of mock experimental runs containing tt 
events, compared to only 0.5% of an ensemble of mock 
experimental runs containing background only. 

We conclude from this sensitivity check that Sleuth 
would not have been able to "discover" tt in the D0 
M^-|-jets data, but that in 30% of an ensemble of mock 
experimental runs Sleuth would have found Vmin la] > 3. 



V. Z+JETS-LIKE FINAL STATES 

In this section we analyze the Z-|-jets-like final states. 
We first describe the data sets and background estimates 
for the dielectron-|-jets channels, and we then discuss the 
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FIG. 6. Scatter plot of where tt Monte Carlo events fall in 
the unit box in the final states W 3j (a) and W ij (c). Al- 
though top quark events appear in the high tails of X^'pt' ^^e 
variable is not particularly discriminating. The locations 
of the data points are shown in (b) and (d). The backgrounds 
are taken to include all standard model processes except top 
quark pair production. 



c 



c 

0) 

.i 

0) 
Q. 
X 
UJ 

o 
o 



o 



35 

30 
25 
20 
15 
10 



Sensitivity check! tl in W 3j(nj) 



w 3j: 

Bkg: 70.5±9.5 
tf: 6.6±2.1 

w 4j: 

Bkg: 9.8±1 .9 
tt: 4.5±1.4 

w 5j: 

Bkg: 1.0±0.3 
tf: 0.73±0.26 

w 6j: 

Bkg: o,i5±o.oe 
tt: o.io±o.04 



-3-2-10 1 2 3 
^ min[G] 

FIG. 7. Histogram of V^in = mm{Vw sj jVw 4] jVw sj , 
Vw&i) for an ensemble of mock experimental runs in which 
the backgrounds include VK+jets and QCD events, and the 
mock samples include (solid) / do not include (dashed) tt in 
addition to the expected background. All experimental runs 
with "Pmin > 3a" are in the rightmost bin. 



dimuon+jets channels. After presenting our results, we 
check the sensitivity of our method to the p resen ce of first 
generation scalar leptoquarks. Appendix B A describes 
signals that might appear in these final states. 



A. Data sets and background estimates 

1. ee2j(nj) 

The ee2j(nj) data set |^l[|, corresponding to an in- 
tegrated luminosity of 123 ± 7 pb~^, is collected with 
triggers requiring the presence of two electromagnetic ob- 
jects. Offline event selection requires two electrons pass- 
ing standard identification criteria with transverse mo- 
menta > 20 GeV and pseudorapidity |?7dct|< 1.1 or 
1.5 <|r7det|< 2.5, and two or more jets with pip > 20 GeV 
and |77dot|< 2.5. At least one electron is required to have 
a matching track in the central tracking detectors and to 
satisfy ionization requirements in the tracking chambers 
and transition radiation detector. For these data the trig- 
ger energy threshold forces a transverse momentum cut of 
20 GeV, rather than the Sleuth-preferred requirement of 
15 GeV. We cut on a likelihood described in Appendix]^ 
in order to correctly identify any events with significant 



missing transverse energy. Electron pairs are combined 
into a Z boson if 82 < rriee < 100 GeV, unless the 
event contains significant (in which case it falls with- 
in eeJpTX, discussed in this section) or a third charged 
lepton (in which case it falls within {£/'y){£/'y){i/ j)X , 
discussed in Sec. 



V]). 



The dominant standard model and instrumental back- 
grounds to this data set are: 

• Drell-Yan + jets production, with Z/j* — > ee; 

• QCD multijets, with two jets faking electrons; and 

• tt pair production with t Wb and with each W 
boson decaying to an electron or to a tau that in 
turn decays to an electron. 

Monte Carlo samples for the Drell-Yan events are gen- 
erated using ISAJET . The Drell-Yan cross section nor- 
malization is fixed by comparing the Monte Carlo events 
with Z + > 2 jets data in the Z boson region. Top quark 
events are generated using herwig at a top quark mass 
of 170 GeV with all dilepton final states included. The 
D0 measured tt production cross section of 5.5±1.8 pb at 
a top quark mass of 173.3 GeV was used |^^. The multi- 
jet background is estimated from a sample of events with 
four or more jets in which the probability for two jets or 
photons to be misidentified as electrons is weighted by 
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Final State 


Z/7*+jets 


QCD fakes 


Total 


Data 


ee2j 


20 ±4 


12.2 ± 1.8 


32 ±4 


32 


ee3j 


2.6 ±0.6 


1.85 ±0.28 


4.5 ±0.6 


4 


edj 


0.40 ± 0.20 


0.24 ± 0.04 


0.64 ±0.20 


3 


ee$T 2j 


3.7 ±0.8 


— 


3.7 ±0.8 


2 


ee$T 3i 


0.45 ±0.13 


— 


0.45 ±0.13 


1 










1 

J. 


Z{-^ ee) 2j 


94± 19 


1.88 ±0.28 


96± 19 


82 


Z{~> ee) 3j 


12.7 ± 2.7 


0.27 ± 0.04 


13.0 ±2.7 


11 


Z{-^ ee) 4j 


1.8 ±0.5 


0.034 ±0.006 


1.8 ±0.5 


1 


Z{-^ ee) 5j 


0.26 ±0.10 


0.0025 ± 0.0009 


0.26 ±0.10 







TABLE IX. Expected backgrounds to the 


ee2j{nj), ee$T2j{nj), and Z{- 


+ ee)2j{nj) final states. 





Final State Z±jets WW ti Total Data 

/i/x2j 0.112 ±0.029 0.25 ±0.13 0.14 ±0.05 0.50 ±0.15 2~ 

fxi^Sj 0.007 ±0.004 0.06 ±0.04 0.065 ± 0.025 0.13 ±0.05 

Z{^ lifi) 2j 2.2 ± 0.4 - 0.050 ± 0.020 2.3 ± 0.4 3 

Z(^^fi)3j 0.24 ±0.05 - 0.018 ±0.009 0.26 ± 0.06 1 



TABLE X. Expected backgrounds to the Z{-~* fifj,) 2j{nj) and fifi2j{nj) final states. 



the number of jets in the event that passed the electron 
Pt and rj requirements. This misidentification probabili- 
ty is calculated from a sample of events with three jets to 
be (3.50±0.35) x 10"'* for an electron with a reconstruct- 
ed track and (1.25 ± 0.13) x 10"'^ for an electron without 
a reconstructed track. The uncertainties in these prob- 
abilities reflect a slight dependence on the jet pt and rj. 
The expected backgrounds for the exclusive final states 
within ee2j(nj) are listed in Table [X. 



2. iJ.iJ.2j{nj) 

The fxfj, 2j{nj) data set ||2^ corresponds to 94 ± 5 pb~* 
of integrated luminosity. The initial sample is composed 
of events passing any of several muon ± jets triggers re- 
quiring a muon with > 5 GeV within |%ot|< 1-7 
and one or more jets with > 8 GeV and |77dot|< 2.5. 
Using standard jet and muon identification criteria, we 
define a final sample containing two or more muons with 
Pt > 20 GeV and |77dct|< 1-7 and at least one muon in 
the central detector (|r/dct|< 1-0), and two or more jets 
with pfj, > 20 GeV and |?7dct|< 2.5. 

We combine a fifi pair into a Z boson if the muon 
momenta can be varied within their resolutions such that 
^fifi ~ Mz and the missing transverse energy becomes 
negligible. More specifically, we combine a muon pair 
into a Z boson if 



X = mm 



/ l/g-l/p^i 
V <5(1/P^i) 



l/b-l/p^2 
<5(1/P^2) 



< 20, 



(4) 



0.7 GeVy^p^/GeV is the error on the missing trans- 
verse energy measured in the calorimeter; niab and fxab 
are the muon pair invariant mass and missing transverse 
energy, computed taking the muons to have scalar mo- 
menta a and &; Mz and Tz are the mass and width of the 
Z boson; and © means addition in quadrature. The cut 
of X < 20 is chosen so that Z{-^ ^fj,) is not the dominant 
background to the fii-j,2j(nj) final states. 

The most significant standard model and instrumental 
backgrounds to this data set are 

• Z + jets production with Z fifi^ 

• WW pair production with each W boson decaying 
to a muon or to a tau that in turn decays to a muon, 
and 

• ti pair production with t — *■ Wb and with each W 
boson decaying to a muon or to a tau that in turn 
decays to a muon. 

A sample oi Z + jets events was generated using VEC- 
BOS, employing HERWIG for parton fragmentation. Back- 
ground due to WW pair production is simulated with 
PYTHIA. Background from tt pair production is simulat- 
ed using HERWIG with a top quark mass of 170 GeV. All 
Monte Carlo samples are processed through a detector 
simulation program based on the GEANT package. 

The expected backgrounds for the exclusive final states 
within iiii2j{nj) are listed in Table |^. The Z{-^ 
fifj,)2j(nj) final states are combined with the Z(—f 
ee) 2j(nj) final states described in S ec. [V A l| to form the 
Z2j{nj) final states treated in Sec. VA3. 



where S{l/p) = 0.18(p - 2)/p'^ 0.003 is the uncertain- 
ty in the reciprocal of the muon momentum; SlfJr) — 
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3. Z2j{nj) 



Combining the results in Tables IX and ^ gives the ex- 
pected backgrounds for the Z2j{nj) final states, shown 
in Table XI . The number of dimuon events in these ta- 



bles is significantly smaller than the number of dielectron 
events due to especially tight identification requirements 
on the muons. 

Z/j* is the dominant background to nearly all final 
states discussed in this section, although other sources 
of background contribute significantly when the dilep- 
ton mass is outside the Z boson mass window. The 
agreement between the total number of events expect- 
ed and the number observed in the data is quite good, 
even for final states with several jets. While any analysis 
of Z-l-jets-like states will need to rely to some degree on 
an accurate Z/7*-|-jets Monte Carlo, having a reliable es- 
timate of the jet distributions in such events is especially 
important when exclusive final states are considered. We 
anticipate that this will become increasingly important 
in the next Tevatron run. Differential agreement between 
data and the expected background may be seen by con- 
sidering a comparison of various kinematic quantities in 
Appendix ^ 



Final State 


Total 


Data 


Z2j 


98± 19 


85 


Z3j 


13.2 ± 2.7 


12 


ZAj 


1.9 ±0.5 


1 


Z5j 


0.26 ±0.10 






TABLE XI. Expected backgrounds to the Z2j{nj) final 
states. 



events in the ee Aj and ee^T 4j final states are provided 
in Appendix 



Data set 



V 



ee2j 
eeSj 
ee4:j 
ee^T 2j 

W2j 
W3j 
Z2j 
Z3j 
ZAj 
Z5j 



0.72 
0.61 
0.04 
0.68 
0.36 
0.06 
0.08 
1.00 
0.52 
0.71 
0.83 
1.00 



TABLE XII. Summary of results on the .Z+jets-Iike final 
states. 



ee 2i 



ee 5j 




B. Results 

The results of applying Sleuth to the Z 2j {nj) and 
M2j{nj) data sets are summarized in Table XII and 
Figs. Ig and ||. Figure || shows the location of the da- 
ta within the unit box for those final states in which the 
two leptons are not combined into a Z boson, while Fig. ^ 
displays the data for those final states in which a Z boson 
has been identified. Large P's are found for most final 
states, as expected. The smallest V's in this class of final 
states are observed in the ee 4j and ee^r 4j final states. 
Although the number of events is small, it is interest- 
ing to compare the number of events observed in the Z 
+ 2, 3, and 4 jet final states (showing good agreement 
with expected backgrounds) with the number of events 
observed in the ee ± 2, 3, and 4 jet and eefJx + 2, 3, and 
4 jet final states. There is a small but statistically in- 
signific ant e xcess in final states with four jets — we find 
in Sec. VII that we expect to find at least one <,0.04 
in the analysis of so many final states. Additionally, one 
of the three ee 4j events has an ee invariant mass barely 
outside the Z boson mass window. The kinematics of the 
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FIG. 8. The positions of the transformed data points in 
the final states ee2j, ee3j, edj, and fifj,2j. The data points 
inside the region chosen by Sleuth are shown as filled circles; 
those outside the region are shown as open circles. 



C. Sensitivity check: leptoquarks 

As a sensitivity check in the Z±jets-like final states 
we consider a scalar, first generation leptoquark |2^ of 
mass rriLQ = 170 GeV, and assume a branching fraction 
to charged leptons oi (3 = 1.0. The cross section for the 
process qq LQLQ with these parameters is 0.54 pb. 
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Sensitivity check! Leptoquarks in ee 2j 



FIG. 9. The positions of the transformed data points in 
the final states Z2j, Z3j, and ZAj. The data points inside 
the region chosen by Sleuth are shown as filled circles; those 
outside the region are shown as open circles. 



The overall efficiency for this type of event is 24±4% , 
including trigger and object requirement efficiencies and 
geometric and kinematic acceptances. If such a lepto- 
quark were to exist, we would expect 11.2 ± 1.5 events of 
signal in the inclusive sample ee2jX, of which 5.9 ± 0.8 
events would fall in the exclusive final state ee2j, on a 
background of 32 ± 4 events. Figure |l^ shows the result 
of Sleuth applied to an ensemble of mock experiments in 
this final state. We see that Sleuth finds V larger than 3.5 
standard deviations in over 80% of these mock samples. 



VI. (£/7)(^/7)(^/7)X 

In this section we analyze the {£/j){£/j){£/j)X final 
states. After describing the data sets and background 
estimates, we provide the results obtained by apply- 
ing Sleuth to these channels. We conclude the section 
with a sensitivity check [X' -> {i/j){£/j){£/-f)X] that is 
more general in nature than those provided for the efj,X, 
VF+jets-like, and Z+jets-like final states above. Exam- 
ples of a few of the many signals that might appear in 
these final states are provided in Appendix B 5 . 
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FIG. 10. Histogram of V for an ensemble of mock experi- 
ments in which the backgrounds include 2/7*-|-jets and QCD 
fakes, and the mock samples include leptoquark pair produc- 
tion (with an assumed leptoquark mass of 170 GeV and P — 1) 
in addition to the expected background. All samples with 
V > 3.5a are in the rightmost bin. Sleuth finds V larger than 
3.5 standard deviations in over 80% of these mock samples. 



A. Data sets and background estimates 

The {£/j){£/j)(£/'-f)X data set corresponds to an in- 
tegrated luminosity of 123 ± 7 pb~^. Global cleanup cuts 
are imposed as above. In this section we strictly ad- 
here to standard particle identification criteria. All ob- 
jects (electrons, photons, muons, and jets) are required 
to have transverse momentum > 15 GeV, to be isolat- 
ed, to be within the fiducial volume of the detector, and 
to be central. For electrons and photons the fiducial re- 
quirement is |77det|< 1-1 or 1.5 <|?7dot|< 2.5; for muons 
it is |?7det|< 1-7. For the case of hadronic jets our cen- 
trality requirement of \r]\< 2.5 is more stringent than the 
fiducial requirement of |77dot| ^4. We require electrons, 
photons, and muons to be separated by at least 0.4 in 
AR = y/{Ariy + (A(/))2. is identified as an object 
if its magnitude is larger than 15 GeV. The selection of 
events is facilitated by use of the database described in 
Ref. 

We make frequent use of the (mis)identification proba- 
bilities determined for these identification criteria, which 



are summarized in Table XIII , 
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Final State Z-y Zj WZ Total Data 

^ 3.3 ±0.7 0.99 ±0.27 - 4.3 ±0.7 3~ 

667 2.1 ±0.4 0.13 ±0.04 - 2.2 ± 0.4 1 

Z-yj 0.80 ±0.30 0.23 ± 0.06 - 1.03 ± 0.31 1 

667j 0.50 ±0.25 0.033 ±0.009 - 0.53 ± 0.25 

ee7^T 0.010 ±0.005 0.024 ± 0.007 0.23 ± 0.10 0.26 ±0.10 1 



TABLE XIV. Expected backgrounds for the 667^ final states. 





6 


7 




6 


0.61 ±0.04 (1^ 


0.28 ±0.03 |30 




7 


0.16 ±0.016^ 


0.73 ± 0.012 |C 




j 


0.00035 ± 0.000035 Q 


0.00125 ±0.00013 


0] 



nu mber of events observed in the data, is given in Ta- 
ble XIV, We find good agreement between the expected 



TABLE XIII. (Mis)identification probabilities. The num- 
ber at (row i, column j) is the probability that the object 
labeling row i will be reconstructed as the object labeling 
column j. 



1. 667X 

The dominant background to eejX is the standard 
model process Z/7*(— » ee)^. We use a matrix ele- 
ment Monte Carlo to estimate this background. The 
pp — > Z/j*{-^ ee)7 cross section, multiplied by our kine- 
matic and geometric acceptance, is 0.50 ±0.05 pb. From 



Table XIII, the probabiUty for two true electrons and 
one true photon to be reconstructed as two electrons and 
one photon is 0.33. From these numbers we estimate the 
expected background from this process into the 667^ fi- 
nal states to be 14.3 ± 2.9 events. Of these, 7.6 ± 1.5 
events have niee < 82 GeV or niee > 100 GeV, and 
82 < mee-y< 100 GcV. Following the prescription in 
Appendix H, such events are placed in the Z final state, 
and are not considered in this section. 

A smaller background in these final states is Z+iets 
production, with the jet faking a photon. From Ref. pq , 
we expect 1100 ± 200 Z{-^ ee)±l jet events in our data; 
the probability that this jet will fake a photon is given 
in Table Kllll. We therefore expect 0.99 ± 0.27 events of 



background in Z7 from this source, and 0.13±0.04 events 
(roughly 10% of the number expected in Z7, determined 
by pythia) in 667. 

The dominant background to the ee^fT final state 
comes from W^(— > ei>)Z{—> ee), in which one of the 
three electrons is reconstructed as a photon. The WZ 
production cross section in the standard model is cal- 
culated to be 2.5 pb jH]; D0's geometric acceptance 
for these events is determined using pyt hia. U sing the 
(mis) identification probabilities in Table XIII , we esti- 
mate the contribution from standard model WZ produc- 
tion to this final state to be 0.23 ± 0.10 events. 

The numbers of expected background events in final 
states with additional jets are obtained by multiplying 
by a factor of 1/5 for each additional jet. The number 
of events expected in each final state, together with the 



background and the numbers of events observed in the 
data. 



2. fJ-fJ-'jX 

The dominant background to the fJ-p-jX final states 
is standard model Z/"f*(—f fJ-fJ-)j. The matrix element 
Monte Carlo used to estimate the backgrounds to 667^ 
is also used for this final state. The normalization 
is determined by multiplying the number of expected 
Z/'y*{-^ ee)7 events by the square of the ratio of efficien- 
cy X acceptance for muons and electrons. For muons, the 
efficiency x acceptance is roughly 0.5 x 0.5; for electrons, 
the number is approximately 0.6 x 0.8. The number of 
expected events in iifij is thus 3.9 ± 0.9. No events are 
seen in this final state. The probability of seeing zero 
events when 3.9 ± 0.9 are expected is 2.8%. 



3. 677X 

The dominant background to e77Ar is the standard 
model process Z/j*{^ ee)7, where on e of th e electrons is 
reconstructed as a pho ton. Fr om Table XIII and the Z(— > 
ee)7 estimate in Sec. VIAl, we determine the number 
of expected events in the e77 final state to be 10.7 ± 
2.1 events. Twelve 677^ events are seen in the data, 
appearing in the final states shown in Table XV. We 
model the e77 backgrounds with the Monte Carlo used 
for the eejX final states above. 

Three of the events in the e77j final state have 
TOe7i72 = 95.8 GeV, TOe7i72 = 85.9 GcV, and rrie-y-^ = 



97.9 GeV, respectively, and are consistent with Zj pro- 
duction with a radiated jet. The invariant masses of the 
objects in the fourth event all lie substantially outside the 
Z boson mass window. Lacking an adequate Z{—f ee)jj 
Monte Carlo, we simply calculate the probability that 
the expected background fluctuates up to or above the 
observed number of events in this final state. The single 



event in the ejj2j final state has rUe 



92.4 GeV; 



this appears to be a boson produced in association 
with two jets. 

One event in this sample contains significant fx in 
addition to one electron and two photons. In this event 
"^•671 = 95.9 GeV, but the missing transverse energy in 
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the event is large, and directly opposite the electron in <f>. 
The transverse mass m^'^ — 71.9 GeV, so this event falls 
in the 11^77 final state. The dominant background to this 
final state is W{-^ ev)Z{-^ ee), in which two electrons 
are reconstructed as photons; the number of such events 
expected in this final state is determined to be 0.11±0.05. 
W{-^ ej^)77 is a slightly smaller but comparable back- 
ground to this final state, which we estimate using a ma- 
trix element Monte Carlo |Q . The total cross section for 
W(-^ ey)"f"f with all three detected objects in the fiducial 
region of the detector and > 15 GeV is determined 
to be 0.77 ± 0.08 fb. The number of W{-^ ev)-i-i events 
in our data is therefore expected to be 0.026 ± 0.010. 
Backgrounds from W^j and W2j, where the jets fake 
photons, are comparable but smaller. This event will be 
combined in the next section with any events containing 
one muon and two photons to form the M^77 final state. 



Final state 


Bkg 


Data 


e77 


10.7 ± 2.1 


6 


> ez/)77 


0.14 ±0.05 


1 


e77J 


2.3 ±0.7 


4 


677 2j 


0.37 ±0.15 


1 



TABLE XV. Population of final states within 677^. 



from this process to be 2.5 ± 0.5 events. The contribu- 
tions from 3j , 7 2j , and 77J are smaller by an order of 
magnitude. 

Two events are seen in the data, both in the final state 
777. One of these events has a three-body invariant mass 
100.4 GeV, consistent with the expectation that 



TO, 



777 



it is truly a Z7 event. The other has a three-body invari- 
ant mass m-yry-y = 153 GeV, but two photons may be cho- 
sen whose two-body invariant mass is rrij^ — 90.3 GeV. 
This event also appears to fit the hypothesis. 



6. eeeX 

The dominant background to the final state eee is 
again Z/7*(— > ee)7, where this time the photon is recon- 
structed as an electron. The cross section quoted above 
for Z/j*{-^ ee)7, fol ded w ith the (mis)identification 
probabilities from Table XII] , predicts 2.6±1.0 events ex- 
pected in the final state eee. One event is seen in the da- 
ta. The eee invariant mass in this event is 87.6 GeV, con- 
sistent with the standard model process Z ee)7, 
where the photon is reconstructed as an electron. 



7. fififiX 



4. fJ.^TjjX 

The dominant backgrounds to the jilpTuX final 
states, like those from the e^TljX final states, come 
from WZ and from a W boson produced in association 
with two photons. The number of expected events from 
WZ is determined as above to be 0.05 ±0.02. The back- 
ground from standard model ^^77 is estimated by multi- 
plying the number of expected W{—^ ev)^^ events above 
by the ratio of efficiency x acceptance for electrons and 
muons. 

Adding the number of events expected from W{-^ 
ei^)77 to the number of events expected from W{-^ 
/ij^)77, we find the total number of expected background 
events in the W^77 final state to be 0.21 ±0.08. No events 
are seen in the muon channel, so the only event in this 
final state is the event in the electron channel described 
above. 



The dominant background to is standard mod- 
el WZ production. We use the WZ production cross 
section above and take our efficiency x acceptance for 
picking up all three muons in the event to be rough- 
ly (0.5 X 0.5)'^ = 0.02. The total number of expected 
background events in /i/x/i from W Z production is thus 
0.020 ± 0.010 events. Zero events are seen in the data. 

The only populated final states within 777X, eeeV, 
and /z^/zV are 777 and eee; these are summarized in 
Table IXVl 



Final state 


Bkg 


Data 


777 


2.5 ±0.5 


2 


eee 


2.6 ± 1.0 


1 



TABLE XVI. Population of final states with three like ob- 
jects. 



5. 777V 



B. Results 



The dominant background to 777 is the standard mod- 
el process Z/^*{-^ ee)j, where both of the electrons are 
reconstructed as photons. Taking th e pro bability of an 
electron faking a photon from Table XIII and using the 
number of Z/7*(— > ee)7 events determined above, we 
find the number of expected events in this final state 



Having estimated the backgrounds to each of these fi- 
nal states, we proceed to apply Sleuth to the data. Large 
V^s are determined for all final states, indicating no hints 
of new physics within {e/-f){e/-f){e/-f)X . Table [XVII 
summarizes the results. We note that all final states 
within {i/j){£/^){£/j)X have been analyzed, including 
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(for example) ee^^lpx and fj,fj,jj2j. All final s tates with- 
in {£/-f){i/j){i/j)X not listed in Table |XVIl| are unpop- 
ulated, and have V — 1.00. 







777 


0.41 


eee 


0.89 


Zj 


0.84 


Zlj 


0.63 


667 


0.88 


ee^IpT 


0.23 


677 


0.66 


e77j 


0.21 


677 2i 


0.30 


W77 


0.18 


TABLE XVII. Summary of results on the 



(lli){ili){lli)X final states. 



effect" in our field to really be a desire for a "3cr effect" 
(one time in one thousand), after a rigorous accounting 
for the number of places that such an effect might appear. 

One of the advantages of Sleuth is that this rigor- 
ous accounting is explicitly performed. The final output 
of Sleuth takes the form of single number, 7^, which is 
"the fraction of hypothetical similar experimental runs 
in which you would see something as interesting as what 
you actually saw in the data." The discussion in the 
preceding paragraph suggests that finding 7^ > Str is as 
improbable (if not more so) as finding a "5cr effect." 

The number of final states that we consider, together 
with the number of background events expected in each, 
defines the mapping between 7'min (the smallest V found 
in any final state) and V . For the final states that we 
have considered in this article, this mapping is shown in 
Fig. [1]. We see that finding V > requires finding 
V > 4.2(7 in some final state. 



C. Sensitivity check: X' {ih){f'h){^h)X 

The backgrounds to the [l li)X final states 
are sufficiently small that a signal present even at the 
level of one or two events can be significant. Due to 
the variety of final states treated in this section and the 
many processes that could produce signals in one or more 
of these final states, our sensitivity check for this section 
is the general process X' — > {I j j j ^^X , rather 
than a specific process such as -pp X2xf ~^ Ui'$T- 
We (pessimistically) take the kinematics of the final state 
particles to be identical to the kinematics of the standard 
model background. In reality the final state objects in 
the signal are expected to have significantly larger mo- 
menta than those in the backgrounds, and the calculated 
V will be correspondingly smaller. With this minimal as- 
sumption about the kinematics of the signal, the details 
of the Sleuth algorithm are irrelevant, and V is given on 
average by the probability that the background fluctu- 
ates up to or above the number of expected background 
events plus the number of expected signal events. 

The quantity V obtained by combining the T^'s cal- 
culated in all final states is a very different measure of 
"significance" than the measure familiar to most high 
energy physicists. The fact that a "significance" of five 
standard deviations is unofficially but generally accepted 
as the threshold for a discovery results from a rough col- 
lective accounting of the number of different places such 
an effect could appear. We can better understand this 
accounting by first noting that five standard deviations 
corresponds to a (one-sided) probability of 3 x 10^''. We 
then estimate that there are at least 5 x 10'^ distinct re- 
gions in the many variable spaces that are considered 
in a multipurpose experiment such as D0 in which one 
could realistically claim to see a signal. A probability of 
1.5 X lO^'^, in turn, corresponds to three standard devia- 
tions. We can therefore understand the desire for a "5fT 
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FIG. 11. Correspondence between "Pmin and V for the final 
states we have considered. 



Let Ny be the smallest integer for which the proba- 
bility that the background in the final state Y fluctuates 
up to or above the expected background b plus Ny is 
< 1.5 X 10~^ (4.2tT). This is the number of events which, 
if observed in y, would correspond to a discovery. This 
number can be related to the most probable cross section 
of the new process into the final state Y through 



Ny 

a^eyC ' 



(5) 



where are the appropriate kinematic and geometric 
acceptance factors for the process d and the D0 detec- 
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tor, cy is the probability that the objects in the true 
final state Y will be correct ly re constructed (which can 
be determined using Table XIII), and £ « 85 pb^^ 



IS 



the effective luminosity of the D0 data after applica- 
tion of global cleanup cuts. The numbers Ny for some 
of the final s tates within {£/'j){£/j){£/j)X are given in 
Table XVIII. (These final states are all unpopulated in 



the D0 data.) Even with our pessimistic assumptions, 
using the Sleuth strategy but setting aside the sophis- 
ticated Sleuth algorithm, we see that a discovery could 
have been made had even a few signal events populated 
one of these channels. 



Final State 


b 


N 




0.059 ± 0.020 


4 


667 2 j 


0.10 ±0.05 


4 


Zl2j 


0.13 ±0.05 


5 


Zl^j 


0.025 ±0.010 


3 


Zj4j 


0.0049 ± 0.0020 


3 




0.10 ±0.05 


4 


e/i/i 


0.040 ± 0.020 


4 




0.020 ±0.010 


3 


W77 


0.21 ±0.08 


5 



TABLE XVIII. The number of signal events A'^ required in 
some of the final states within {£/'y){£/y){£/'y)X in order to 
find V > 3a (see the discussion in the text). This number is 
pessimistic, as it assumes that the signal is distributed iden- 
tically to the backgrounds in the variables of interest. Most 
tenable models predict events containing final state objects 
that are significantly more energetic than the backgrounds, 
and in this case A'' decreases accordingly. 



VII. SUMMARY 



Table KIX summarizes the values of V obtained for 
all populated final states analyzed in this article. Taking 
into account the many final states (both populated and 
unpopulated) that have been considered in this analy- 
sis, we find V — 0.89 (— 1.23cr). Figure ^ shows a his- 
togram of the V^s computed for the populated final states 
analyzed in this article, together with the distribution 
expected from a simulation of many mock experimental 
runs. Good agreement is observed. 

Although no statistically significant indications of new 
physics are observed in this analysis, some final states 
appear to hold greater promise than others. The smallest 
■p's (0.04 and 0.06) are found in the final states eeAj 
and ee^T 4j. The kinematics of the events in these final 
states are provided in Appendix 

It is very difficult to quantify the sensitivity of Sleuth 
to arbitrary new physics, since the sensitivity necessar- 
ily depends on the characteristics of that new physics. 
We have provided examples of Sleuth's performance on 
"typical," particular signatures. This function is served 
by the sensitivity checks provided at the end of each of 



Data set 






V 














0.14 


(±1.08cr) 






0.45 


(±0.13cr) 


efj,^T 2 j 




0.31 


(±0.50cr) 




VF±jets-like 


0.71 


(-0.55cr) 


W2j 




0.29 


(±0.55cr) 


W3j 




0.23 


(±0.74cr) 


WAj 




0.53 


(-O.OScr) 


W5j 




0.81 


(-0.88cr) 


W6j 




0.22 


(±0.77(j) 


e$T 2j 




0.76 


(-0.71cr) 


eJ^T 3j 




0.17 


(±0.95cr) 


6^T 4j 


Z±jets-like 


0.13 


(±1.13cr) 


Z2j 




0.52 


(-0.05cr) 


Z3j 




0.71 


(-0.55cr) 


Z4j 




0.83 


(-0.95cr) 


ee2j 




0.72 


(-0.58a) 


66 3j 




0.61 


(-0.28cr) 


66 4j 




0.04 


(±1.75cr) 


66_^T 2j 




0.68 


(-0.47cr) 


eefJr 3j 




0.36 


(±0.36cr) 


66_^T 4j 




0.06 


(±1.55cr) 


fj,fi2j 


(£h){£h){£h)X 


0.08 


(±1.41cr) 


eee 




0.89 


(-1.23cr) 


Zj 




0.84 


(-0.99cr) 


ZlJ 




0.63 


(-0.33cr) 


667 




0.88 


(-1.17cr) 


ee^IpT 




0.23 


(±0.74cr) 


677 




0.66 


(-0.41cr) 


6777 




0.21 


(±0.81cr) 


677 2j 




0.30 


(±0.52cr) 


Wjj 




0.18 


(±0.92cr) 


777 




0.41 


(±0.23cr) 


V 




0.89 


(-1.230-) 



TABLE XIX. Summary of results for populated final 
states. The most interesting final state is found to be eeAj, 
with V — 0.04. Upon taking into account the many final 
states we have considered using the curve in Fig. |l^, we find 
V = 0.89. The values of V obtained in these final states are 
histogrammed in Fig. |l^, and compared to the distribution 
we expect from an ensemble of mock experimental runs. No 
evidence for new high pr physics is observed in these data. 



Sees. ^-|vi| In the analysis of the efiX data in Rcf. [||, 
our signal was first WW and tt together, and then only 
tt. This was a difficult signal to find, for although both 
WW and tt cluster in the upper right-hand corner of 
the unit box, as desired, we expect only 3.9 WW events 
in e^fJr (with a background of 45.6 events), and 1.8 ti 
events in e/i^r 2j (with a background of 3.4 events). We 
were able to consistently find indications of the presence 
of WW and tt in an ensemble of mock experiments, but 
we would not have been sufficiently sensitive to claim a 
discovery. 
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VIII. CONCLUSIONS 




-4 -5 -2 -1 

FIG. 12. Histogram of the P's computed for the populated 
final states considered in this article. The distribution agrees 
well with the expectation. 



We have applied the Sleuth algorithm to search for 
new high pT physics in data spanning over thirty-two 
exclusive final states collected by the D0 experiment 
during Run I of the Fermilab Tevatron. A quasi-model- 
independent, systematic search of these data has pro- 
duced no evidence of physics beyond the standard model. 
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In the W^-|-jets-like final states we again chose tt for 
our sensitivity check. This was both a natural sequel 
to the sensitivity check in efiX and a test of Sleuth's 
performance when the signal populates the high tails of 
only a subset of the variables considered. We find Vmin > 
3(7 in 30% of an ensemble of mock experimental runs 
containing tt events on the final states W3j, W 4j, W 5j, 
and W6j, compared with only 0.5% of an ensemble of 
mock experimental runs containing background only. 

In the Z-t-jets-like final states we considered a lepto- 
quark signal. This is in many ways an ideal signature — 
a relatively large number of events (about six) are pre- 
dicted, and the signal appears in the high tails of both 
variables under consideration. Sleuth finds V > 3.5cr in 
over 80% of the mock experiments performed. 

Finally, in the final states {£/'y){£/'y){£/'y)X we intro- 
duced the mapping between Vmin and V and briefly dis- 
cussed its interpretation. The generic sensitivity check 
we considered [X' {l/j){£/-f){£/-f)X] demonstrates 
the advantages of considering exclusive final states. 
While the other sensitivity checks rely heavily upon the 
Sleuth algorithm, this check shows that a careful and 
systematic definition of final states by itself can lead to 
a discovery with only a few events. 



APPENDIX A: DEFINITIONS OF FINAL STATES 

This appendix reviews the definitions of final states 
provided in Ref. |3j . The specification of the final states 
is based on the notions of exclusive channels and stan- 
dard particle identification. We partition the data into 
exclusive final states because the presence of an extra 
object (electron, photon, muon, . . . ) in an event often 
qualitatively changes the probable interpretation of the 
event and the variables that naturally characterize the fi- 
nal state, and because using inclusive final states can lead 
to ambiguities when different channels are combined. 

We attempt to label these exclusive final states as com- 
pletely as possible while maintaining a high degree of 
confidence in the label. We consider a final state to be 
described by the number of isolated electrons, muons, 
photons, and jets observed in the event, and whether 
there is a significant imbalance in transverse momentum. 
We treat as an object in its own right, which must 
pass certain quality criteria. In Run I D0 was unable to 
efficiently differentiate among jets arising from b quarks, 
c quarks, light quarks, and hadronic tau decays. We con- 
sider final states that are related through global charge 
conjugation to be equivalent in pp or e^e~ (but not pp) 
collisions. Thus in principle e+e^7 is a different final 
state than e^e^7, but e~''e^7 and e~e~7 together make 
up a single final state. D0 lacked a central magnetic field 
in Run I, so we choose not to distinguish between e+/e^ 
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or fj,~^ In events containing two same-flavor leptons, 
we assume that they are of opposite charge. 

We combine an e+e" pair into a Z boson if their in- 
variant mass me+e- fahs within a Z boson mass window 
(82 < TTie+e- < 100 GcV) and the event contains nei- 
ther significant nor a third charged lepton. A fJ-'^ fi~ 
pair is combined into a Z boson if the event can be fit 
to the hypothesis that the two muons are decay products 
of a Z boson and that the in the event is neghgible, 
and if the event contains no additional charged lepton. 
If the event contains exactly one photon in addition to a 
i~^£~ pair, and contains neither significant nor a third 
charged lepton, and if nig+g- does not fall within the Z 
boson mass window, but mi+i-^ does, then the £'^£~j 
triplet becomes a Z boson. An electron and become 
a W boson if the transverse mass m^p, is within a W 

boson mass window (30 < rri^^^ < 110 GeV) and the 
event contains no second charged lepton. A muon and 
in an event with no second charged lepton are always 
combined into a W boson; due to our more modest muon 
momentum resolution, no mass window is imposed. Be- 
cause the W boson mass window is so much wider than 
the Z boson mass window, no attempt is made to identify 
radiative W boson decays. 



APPENDIX B: EXAMPLES OF SIGNALS THAT 
MIGHT APPEAR 

In this section we provide a few examples of signals 
that might have been discovered in the course of this 
analysis. This discussion is provided to give the reader a 
taste of the many processes that might appear in the final 
states we have analyzed, and is by no means intended to 
be complete. The possibility that the correct answer is 
"none of the following" is one of the strongest motivations 
for pursuing a quasi-model-independent search. 



1. efj,X 



In supersymmetric models (denoting the supersym- 
metric particles as in Ref. [Q), the process qq Z/j* ~* 
Xi^Xi ~* e/i^j^x^Xi "^^-n produce events appearing in the 
e^fx final state. More generally, any process involving 
the production of two charginos has the potential for pro- 
ducing a final state containing an electron, a muon, and 
^T- This final state may also be reached through the lep- 
tonic decays of two taus, obtained (for example) from the 
production of two f particles that each decay to txI, or 
from the production of a heavy Z-like object that couples 
strongly to the third generation. An anomalous correc- 
tion to the standard model WW^ vertex or anomalies 
involving the top quark could also appear in these final 
states. 



2. Final states already considered 

A sampling of the types of new physics that might 
appear in a few of the final states described in Sec. [II A 
is provided here. 

2j. The dijet final state could contain hints of a mas- 
sive object (such as an additional neutral gauge boson) 
produced through qq annihilation and decaying back in- 
to qq. It could also contain indications that quarks 
are in fact composite objects, interacting through terms 
in an effective Lagrangian of the form -^qqq'q', where 
A ^ 1 TeV is a compositeness scale and c is a constant of 
order unity. 

eflT- Models containing symmetry groups larger than 
the SU(3)cxSU(2)i;,xU(l)Y group of the standard mod- 
el often contain an additional SU(2) group, suggesting 
the existence of a heavy W^-like gauge boson (W) that 
would decay into the efx final state, with the transverse 
mass of the electron and neutrino greater than that ex- 
pected for the standard model W. Production of ID de- 
caying to ixii^Xi could also produce events in this final 
state, as could production of xfx2 decaying to ivxii^^Xi- 

ee. If both quarks and leptons are composite ob- 
jects, there will be four-fermion contact terms of the form 
-^qqt'^t~ in addition to the -^qqq'q' terms postulated 
in the discussion of the 2j final state above. Such an 
interaction would produce events with large transverse 
momentum, opposite-sign leptons, and should appear in 
the ee and /i/i final states. Some models that employ a 
strong dynamics to break electroweak symmetry predict 
the existence of composite "techni-" particles, such as the 
LOT, Pt, and ttt, that are analogous to the composite w, p, 
and TT mesons that arise from confinement in QCD. The 
technirho (px) and techniomega (ujt), if produced, will 
decay into an £'^£~ pair if their preferred decay mode to 
technipions (ttt) is kinematically forbidden. Such events 
will appear as a bump in the tail of the ee invariant mass 
distribution and as an excess in the tail of the electron pT 
distribution. Models containing symmetry groups larg- 
er than that of the standard model typically contain a 
heavy neutral boson (generically called a Z') in addition 
to the W boson described above. If this Z' boson cou- 
ples to leptons, the process qq ^ Z' ^ ££ could produce 
a signature similar to that expected from the decay of a 
Pt or LOT- 



3. VF+jets-like final states 

A variety of new signals have been predicted that 
would manifest themselves in the VF-|-jets-like final states 
— those final states containing events with a single lep- 
ton, missing transverse energy, and zero or more jets. A 
plethora of supersymmetric signatures could appear in 
these states. A chargino and neutralino, produced from 
qq through an s-channel W boson, can proceed to decay 
as Xi^ ^^Xi and X2 ~^ QQXii leaving an event that 
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will be partitioned into either the elpT 2j or W 2j final 
state. Pair production of top squarks, with i — + bx^ 
and subsequent decays of the charginos to ei^Xi ^^id 
QQ'Xii ■will produce events likely to fall into the e^r -^j or 
W 4j final states. Depending upon the particular model, 
even gluino decays can give rise to leptons. Events with 
gluinos that are pair-produced and decay, one into qq'xf 
and the other into qqx\, can also find themselves in the 
eJpT^j or VF4j final state. Other possible decays of the 
supersymmetric spectrum allow many more signals that 
might populate these final states. 

The decay of a p^, produced by qq annihilation, can 
produce a boson and a tt^ , which in turn may decay 
to hb or gg. Such an event should appear in the high tails 
of the and plj, distributions in our analysis of the 
W 2j final state if the technipion is sufficiently massive. 
The same final state may also be reached by the pro- 
cess gg ^ — > W~7T^ — > £~i^cb. A neutral color-octet 
technirho (pxs) produced by gg annihilation can decay to 
two technipions carrying both color and lepton quantum 
numbers (ttlq), each of which in turn decays preferen- 
tially into a massive quark and a massive lepton. If the 
technipion is heavier than the top quark then the decay 
ttlq — > tr or ti>r is kinematically allowed. Appropriate 
decays of the W bosons from the two top quarks leave the 
event containing one high transverse momentum lepton, 
substantial ]^t, and several energetic jets. 

The standard model contains three generations of 
quarks and leptons, but there appears to be no fun- 
damental reason that Nature should choose to stop at 
three. A massive charge —1/3 fourth-generation quark 
(6'), which could be pair-produced at the Tevatron, would 
be apt to decay weakly into a W boson and a top quark. 
Events in which one of the four W bosons then decays lep- 
tonically will result in a final state containing one lepton, 
substantial missing transverse energy, and many jets. 

Leptoquarks, a consequence of many theories that at- 
tempt to explain the peculiar symmetry between quarks 
and leptons in the standard model, could also be pair- 
produced at the Tevatron. If their branching ratio to 
charged leptons f3 = 0.5 then the pair will decay to ivqq 
50% of the time, resulting in events that will be classified 
either as e^r^j or W2j. 

Models invoking two Higgs doublets predict a charged 
Higgs that may appear in occasional decays of the top 
quark. In such models a top quark pair, produced by gg 
or gg annihilation, can decay into H~^bW~b. Depend- 
ing upon the mass of the charged Higgs particle, it may 
decay into W'^bb, cs, or t^v. Appropriate decay of the 
W boson(s) in the event will result in the event populat- 
ing one of the W2j{nj) final states. Other predictions 
abound. 



4. Z+jets-like final states 

Just as in the Il^-|-jets-like final states, there are a host 
of theoretical possibilities for new physics in the Z-|-jets- 
like final states. Although some of these processes involve 
the production of two same-flavor, opposite-sign leptons 
via the production of a standard model Z boson, many 
others involve particles that decay to leptons of different 
flavor, or with the same charge. These different possi- 
bilities typically are partitioned into different final states 
according to our prescription: events that contain leptons 
of d iffere nt flavor (those within ejiX) are considered in 
Sec. II C; events containing leptons of similar charge (e.g.. 



an e"^e+ 2j event) would in principle be partitioned into 
different flnal states than events containing leptons of op- 
posite charge (e.g., an e+e" 2j event) if D0 distinguished 
electron charge; and events in which the leptons have an 
invariant mass consistent with the hypothesis that they 
are the decay products of a Z boson are partitioned into 
different flnal states than those with a dilepton invariant 
mass outside the Z boson mass window. 

Models containing supersymmetry and imposing con- 
servation of i?-parity predict signatures containing sub- 
stantial missing transverse energy. Such events might 
therefore populate the eeflr 2j(nj) or ^i^fr 2j(nj) chan- 
nels. Final state leptons may be obtained in supersym- 
metric models from the decays of neutralinos (which can 
produce two same-flavor, oppositely-charged leptons), or 
charginos or sleptons (which decay into a single charged 
lepton and missing transverse energy). The process 
gg' W* ~* Xi X21 'with subsequent decay of the 
chargino to gg'xi and the ncutralino to £^i~Xij results 
in an event with two same-flavor, opposite-sign leptons, 
two jets, and missing transverse energy, and would ap- 
pear in our eefx 2j or fj,fj, 2j flnal states. Events in which 
gluinos are pair-produced and decay via g — s- qq'xt will 
appear in the eefx^^j and fJ-fJ-'ij flnal states when the 
gaugino decays to ivxi- Pair production of scalar top 
quarks {qq/gg 9 ^ it*) that decay via i —>■ bxf 
and xf —> ii^Xi again produce events that populate 
the eelpT 2j and /i/z 2j flnal states, in addition to the 
e^fx 2j flnal states already considered. If i?-parity is vi- 
olated, then supersymmetric signals could populate flnal 
states without missing transverse energy. Pair produc- 
tion of gluinos decaying to ccl could produce events that 
land in the ee 4j flnal state if the i?-parity- violating decay 
cl — > e^d is allowed. 

Color-octet models predict the existence of a color- 
octet technirho, which can decay to TTLgT^LQ. These tech- 
nipions decay preferentially to massive particles, like the 
color-singlet ttt, but their decay products will carry both 
color and lepton quantum numbers. Events in which each 
TTLQ decays to a & quark and a r lepton will populate 
eelpT 2j and /i/i 2j flnal states, among others. Lepto- 
quarks motivated by grand uniflcd theories could be pair- 
produced at the Tevatron via gg — > Z/"f* LQLQ, and 
might populate the flnal states ee2j and ^^i2j. Again, 
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other examples abound. 



5. {ihWhW/^)x 

There are few standard model processes that produce 
events in which the sum of the numbers of electrons, mu- 
ons, and photons is > 3. The {ih){ih){lh)X final 
states are therefore quite clean, and the presence of even 
a few events in any of these states could provide a strong 
indication of new physics. 

Supersymmetric models predict a variety of possible 
signatures in these states. Those models in which R- 
parity is conserved produce events with missing trans- 
verse energy in addition to three {l/^) objects. Models 
in which the lightest neutralino (x?) is the lightest su- 
persymmetric particle (LSP) usually produce final states 
without photons. This case occurs for many models in 
which the supersymmetry is broken in a hidden sector 
and communicated to the visible sector through gravi- 
tational forces (gravity-mediated supersymmetry break- 
ing). Models in which the gravitino (G) is the LSP often 
produce final states with photons from the decay of the 
next-to-lightest supersymmetric particle via (for exam- 
ple) Xi ~* iG- This case, in turn, obtains for many 
models in which the breaking of the supersymmetry is 
mediated by gauge fields (gauge-mediated supersymme- 
try breaking). For example, the production of a chargino 
and neutralino through qq annihilation into a virtual W 
boson can produce events in these final states through 
the decays xf f-^Xi ^^^d X2 ^^Xi if the lightest neu- 
tralino is the LSP, or through the decays Xi ~^ ci'Xij 
X2 ~^ QQXij and x5 ~' 7^ if the gravitino is the LSP. 

Charginos can be pair-produced in the reaction qq 
Z/7* — > xfxf ■ If they decay to eux2, and if X2 in turn 
decays to 7X1, these events will populate the final state 
ee77^T- The production of slepton pairs can also result 
in events falling into the final state eejjfx, since a typi- 
cal decay of a selectron in a model with gravity-mediated 
supersymmetry breaking is e ^ 6X21 with 

X°2 ^ 7X5 ■ If 

a pair of sufficiently massive sleptons are produced, each 
can decay into the corresponding standard model lepton 
and the second- lightest neutralino (X2)) which in turn 
could decay into Uxi- A similar production of ii> can 
easily lead to a final state with one fewer charged lep- 
ton, through the decay chain i> —f £x^ , and xf — > ^i^Xi- 
The standard model backgrounds to such events, contain- 
ing five or more charged leptons and substantial missing 
transverse energy, are vanishingly small. Events with 
four charged leptons and substantial fix could result from 
the decay of a X2X2 pair, in which each x" decays to iixi- 
Even pair production of gluinos, each decaying to qqX2^ 
with one neutralino decaying to eex? and the other to 
7X?, could produce events in these final states. With 
this particular decay, such events would appear in the 
final state eej 2j. 

If leptons exist in excited states several hundred GeV 



above their ground state, just as hadrons exist in excited 
states at energy scales a thousand times smaller, they 
could be produced in the process qq — > Z/7* £*£* or 
qq' — > W* t*v*. The excited leptons can decay by 
emitting a photon, so that f — > £7 and v* — + v^. Such 
events would populate the ^£77 and ifxll final states. 
If the technirho exists and is sufficiently massive, it can 
decay to WZ . Roughly one time in fifty both the W and 
Z bosons will decay to leptons, producing a £'^£~t'$T 
event. More generally, any process producing anomalous 
triboson couplings will affect t he (£/ 'y)(£/'y)(£/'y)X final 
states, and (as we show in Sec. VI C[ ) our method is likely 
to be sensitive to such a signal. 



APPENDIX C: COMPARISON OF 
DISTRIBUTIONS 

In this appendix we show kinematic distributions of 
the data and expected backgrounds for the most heavi- 
ly populated final states that we have considered. Fig- 
ures p^|-p^ show good agreement between data and the 
expected background in a number of distributions for the 
heavily populated VF-|-jets-like final states W2j, W3j, 
and W 4j . Figures |l^ and |lj serve the same function for 
the final states Z 2j and Z3j. 
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FIG. 13. Comparison of background to data for W2j. 
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APPENDIX D: SIGNIFICANCE 

We determine the significance of any missing trans- 
verse energy in an event in the Z+jets-hke final states 
by computing a probabihty density pifr)- This is a true 



probabihty density in the sense that, for a given event, 
the probability that the actual missing transverse ener- 
gy in that event is between and ^t + S^t is given 
by p{$t)5$t- This density is computed with a Monte 
Carlo calculation. For each data event we generate an 
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ensemble of events similar to the original but with the 
energies of the objects smeared according to their res- 
olutions. Jets are smeared with a Gaussian of width 
tJ = 80%VE, and electrons are smeared with a Gaussian 
of width cr ~ 20%y/E (a slight inflation of the measured 
resolution of lb%VE), where E is the energy of the ob- 
ject in GeV. The component of the missing transverse 
energy along the direction of the original is re- 
calculated for each smeared event, and the values that 
are obtained are histogrammed. The histogram is then 
smoothed, and the likelihood 



is calculated. 



P{$T, = 0) 

Studies have shown 



(Dl) 



that a cut of 



logj^o -C^T > 3 does an excellent job of retaining events 
with true while rejecting QCD background. 



APPENDIX E: KINEMATICS OF INTERESTING 
EVENTS 



Table pCX| provides information about the events in the 
most interesting final states seen in the course of this 
analysis. Invariant masses of objects in these events are 
given in Table XXI, 
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TABLE XX. Kinematic properties of the most interesting 


events seen in this analysis. 
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TABLE XXI. Invariant masses (in units of GeV) of objects 
in the most interesting events seen in this analysis. 
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